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Abstract 

We propose an effective exponent ruling the algebraic decay of 
the average quantum return probability for discrete Schrodinger op- 
erators. We compute it for some non-periodic substitution potentials 
with different degrees of randomness, and do not find a complete qual- 
itative agreement with the spectral type of the substitution sequences 
themselves, i.e., more random the sequence smaller such exponent. 
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Anomalous transport in non-periodic structures is due to intricate quan- 
tum interferences which may also lead to localization of wave functions. An- 
other possibility is ballistic motion, mainly related to periodic structures. 
Here we consider transport properties in nearest-neighbours tight-binding 
models in 7L whose general Hamiltonian H is given by 

{Hij)) n = ip n+1 + Ipn-! + \V n ip n , (1) 

with A > and potentials V = (V^)„ g 2 generated by some non-periodic 
substitution sequences. 

Among the characterizations of (de) localization and transport we single 
out the (average) moments of the "position" operator 

I rT 00 

m a (T) = - dt J2 \n-n \ a \Mt)\ 2 , a^O, (2) 

1 J ° n=-oo 

and the (average) return probability 

C(T) = ±J o T dt\^ no (t)\ 2 - (3) 

In relations (§) and (§) it is implicitly assumed that the initial condition is 
i>n — ^n,n - Both quantities m a (T) and C(T) have strong physical appeal and 
in some cases are attainable to theoretical and numerical investigations. We 
notice that the return probability was one of the first quantities considered 
in the seminal paper by Anderson on localization in disordered structures 
It has been found that for large T 0, |, |, |, |, § 

m a (T) ~ T Q/3(a) and C{T) ~ T~ A . (4) 

Localization should be characterized by vanishing exponents f3 and A, bal- 
listic motion by (3(a) = A = 1, while anomalous transport by < /3, A < 1. 
Notice that (3(2) is related to the direct conductivity via the anomalous Drude 
formula g |. 

In this work we consider potentials V in ([]]) generated by some substitu- 
tion sequences and compute the exponent for the decay of the return proba- 
bility as a function of the degree of randomness of those sequences. As it will 
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be seen later, in most cases A can not be obtained directly from numerical in- 
tegration of the time-dependent Schrodinger equation; namely, the standard 
fitting procedure given by equation @ works only for the well-investigated 
case of Fibonacci potentials (to be defined below). We then propose an al- 
ternative approach, based on the energy spectral decomposition of the initial 
quantum state, which not only is able to retrieve the known Fibonacci results 
but permits us to exhibit an effective exponent also for other non-periodic 
substitution potentials such as Thue-Morse, Rudin-Shapiro, paperfolding and 
period doubling (see below for their precise descriptions). The main conclu- 
sion will be that there is no perfect correspondence of the exponent ruling the 
algebraic decay of the return probability and the degree of randomness of the 
own substitution sequences, as one expects based on general considerations 
(i.e., more random the sequence smaller such exponent). Now we proceed to 
the details of the points just outlined, including a justification for the choice 
of A as our exponent of interest. 

The RAGE theorem and Wiener lemma give direct physical meaning to 
the standard types of spectra, i.e., point and continuous (absolutely and sin- 
gular) in the sense that one of their corollaries is that for continuous spectra 
(3(a) > and A > 0, while for point spectra A = 0; it is left room for 
(3 > even for point spectra due to the tail of eigenfunctions and domain 



intricateness |10|, IT]. The determination of such exponents is a quantitative 



step from RAGE and Wiener lemma which, by its turn, is still related to 
deep spectral quantities, i.e., generalized dimensions of the (positive) spec- 
tral measures associated to the initial state ip. It has been rigorously proven 
that (3(a) is bounded from below by the information dimension D\ |2|, [3| (all 
dimensions are related to the corresponding spectral measure) and also con- 
jectured that (3(2) ~ D Q (D denotes the fractal dimension of the spectrum). 
In an interesting paper Guarneri and Mantica |T2[ have presented examples 
of homogeneous fractal spectral measures, i.e., with generalized dimensions 
D q = D for any q, for which (3(a) is not constant and no simple exact re- 
lation seems to hold between the thermodynamics of the spectrum and the 
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exponent /3, so that "multiscaling does not require multifractality" P^j ; this 
was called quantum intermitency in |12, 13 . See |14[ for some recent results 
on a particular class of systems and other references. 

On the other hand, it was rigorously proven |5], || that the exponent A 
ruling the algebraic decay of the return probability equals the correlation 
dimension D 2 . We note that such relation supposes the limit defining the 
dimension D 2 does exist. See also M where this relation was first proposed 
in the context of anomalous diffusion. Therefore we have selected the return 
probability, its corresponding exponent A and correlation dimension D 2 , as 
the main tools for analyzing our systems. 

Relevant examples of anomalous diffusion are generated by almost-periodic 
potentials V; an important class of such potentials is given by V induced by 
non-periodic (primitive) substitution sequences [[15], jlEfl . These sequences 
form a convenient laboratory for the study of anomalous transport since in 
all rigorously analyzed cases they generate singular continuous spectra for 
the tight-binding model (|1|), although the own spectral types of substitution 
sequences are not equal; for example, Fibonacci (Fee), paperfolding (PF) 
and period doubling (PD) substitution sequences have point autocorrelation 
measures, Thue-Morse (TM) has singular continuous autocorrelation mea- 
sure, and the autocorrelation measure of Rudin-Shapiro (RS) substitution 
sequence is absolutely continuous. Although all these sequences are almost 
periodic, their spectral properties characterize them qualitatively from "or- 
dered to random" since periodic and quasiperiodic sequences have pure point 
autocorrelation measures (as Fee, PF and PD do), whereas independent 
random sequences have absolutely continuous autocorrelation measures (as 
RS does). The TM sequence lies in an intermediate place. It is also worth 
noting that all these sequences give rise to strictly ergodic dynamical systems 
with zero topological and generalized entropy JT^ |T7|j . 

Due to different degrees of randomness of the substitution sequences, it 
was expected differences in the spectral properties of the corresponding tight- 
binding Hamiltonians ([J), but as already commented above all rigorously 
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studied cases have presented singular continuous spectrum |H| [H], [21], |^T|. 



22J, m (the RS case is an important open problem 0, g|, 05 the PF 



case is also open). Our main goal in this Letter is to investigate whether 
the spectral character of the sequence generating the potential is responsible 
for different physics through details of the return probability behaviour. To 
this end we consider potentials V in ([3]) generated by the five non-periodic 
sequences Fee, PF, PD, TM and RS. A can be computed either from 
numerical integration of the time-dependent Schrodinger equation, and then 
fitting a straight line on logC(t) x logt, or directly computing D2 from its 
definition (which also involves a straight line fitting — see below). However, 
as already anticipated, we have faced problems in linear fittings in both 
procedures (except for the well-investigated case of Fee potentials) and we 
propose a pragmatic approach to get such exponents which is able to recover 
the Fee known results. 

Now we present the rules describing the sequences we use to generate V. 
Fee, PD and TM sequences are constructed with an alphabet of two letters 
{a, b} through the substitutions 



a — > ab, b — > a (Fee), a — > ab, b —> ba (TM), 

a — > ab, b — > aa (PD). 

Beginning with a and applying successively the substitution rules, non-periodic 
sequences are obtained; e.g., the Thue- Morse sequence is given by 

abbabaabbaababba ■ ■ ■ 

The RS and PF sequences can be obtained with an alphabet of four letters 
{a', b', c', d'}, the substitutions 

a'^a'b', b'^a'c', c'^d'b', d' -»• d'e' (RS), 

a'^a'b', b'^c'b', c'^a'd', d! -»• c'd' (PF), 
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and then the identifications a', b' — > a and c', d' — > b in both cases; the first 
elements of the RS sequence are 

aaabaabaaaabbb ■ ■ ■ 

We then use these substitution sequences to define our potentials V; we 
take V n = if the n-th letter of the sequence is a and V n — 1 in case it 
is b. There are standard ways to extend the potential for negative values 



of n [fit], |22|], but we avoid such issue by taking a finite sample of N sites, 
with n > 0, and using the initial wavefunction ip n = 5N/2,n - 111 this way we 
construct the almost-periodic substitution potentials AV^ and investigate A 
as function of the degree of randomness of the potential and its intensity A. 

It is known that Fee, PD and TM generate potentials whose spectra 
of (U) are singular continuous for all A 7^ 0. The case of RS has been 



numerically investigated in [T2^, [2(J indicating point spectrum for A > 2 and 
mixed spectrum, i.e., point and singular continuous, for < A < 2 (notice 
we use a scale for the potential values which is different from [Pq| ). For 
Hamiltonian (Jl|) with PF potential it is only known that its spectrum has no 



absolutely continuous component, since it is primitive p2| ; from a rigorous 
point of view the lack/presence of eigenvalues in this case is also an open 
question. 

The case of Fee Hamiltonian has also been considered in |4[] and a good 
agreement between the value of A from numerical integration of Schrodinger 
equation and D 2 was found. Let's recall the definition of D 2 associated to a 
spectral measure /1 and how it is usually estimated 0, ||. For e > let B £ (x) 
denote the open ball of centre x and radius e/2 and set 



7(e) = / ii{B £ {x))dn{x)\ 
the correlation dimension of fi is given by the limit 



B 2 = llm !^. (5) 

£^0 log£ 
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If this limit does not exist one defines D 2 and D 2 via limsup and liminf, 
respectively. The latter remark is important here since we have found numer- 
ical indications that the limit in the definition (|j) does not exist for general 
substitution potentials. In numerical practice we have a finite basis approxi- 
mation for (P whose spectrum is composed of eigenvalues Xk ] then we divide 
the energy range into boxes Bj of length e, approximate 

7(^7*(^) = E( E M 2 ) 2 (6) 

and get D 2 from the linear fitting of log7*(£) x loge. is the projection of 
the initial wavefunction ip on the eigenvector with eigenvalue Xk- We have 
used these procedures to recover D 2 and A found in || for the Fee case as 
illustrated in figure 1. However such techniques do not work for substitution 
potentials distinct from FCC (as far as we have checked), since no clear region 
with linear behavior is found in the plots log7*(£) x loge and logC(t) x logt, 
as exemplified in figure 2 for the PD potential with A = 1.6. We suspect this 
behaviour is an indication that the limits defining the scale exponents D 2 
and A are not well defined in such situations; then we propose a pragmatic 
approach to extract effective exponents D 2 by selecting a particular value e* 
of e. Before presenting our approach we stress we have also tried to get 
well-defined exponents by site averaging on samples beginning at locations 
0, 1 x 10 4 , 2 x 10 4 , • • • , 5 x 10 4 , but quite similar behaviours were found. 

We begin our argument with the remark that if e is smaller than the least 
eigenvalue spacing (we just ignore the possibility of degenerate eigenvalues 
in this argument) then 

7*(^) = EKI 4 

k 

which resembles the so-called inverse participation ratio (which does not 
depend on e); this also gives a physical interpretation for D 2 . As a naive 
first guess for an effective exponent one could try to use J2k l°fc| 4 instead of 
7*(e), but the exact value of e to be used in an approximation to fl5|) is not 
clear at all. The theoretical determination of D 2 involves the limit e — > 0; 
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finite basis approximations preclude this limit and also too small values of 
e are meaningless, despite the inverse participation ratio interpretation. For 
sufficient small values of e we have 

-^L«l- (7) 

V \n, I 4 1 ' 

l^k \ a k\ 

we suggest to pick e* as the smallest value of e such that the l.h.s. of (0) 
considerably deviates from 1, so still keeping track of 7*(e) and also the 
inverse participation ratio interpretation in operation. Then we estimate the 
effective D-i as D\ given by 

2 loge* v ' 

Let's be more precise on how we have picked up e* in practice. By using 
double precision (16 digits) in our code, we adopted that after diagonalization 
we can numerically resolve the spectral quantities with 8 digits, i.e., half of 
the number of digits of the code precision, so that e* is given by the smallest 
value of e such that 7*(e) — J2k \ a k\ 4 > 1CT 8 or, equivalently, the smallest e 
such that 

7*0) , 10- 8 
- 1 > 



EfcKI 4 EfcKI 4 

We remark that in most cases e* can also be obtained directly from visual 
inspection, as in figure 2b, and the precise value 10~ 8 is not so relevant since 
in general j*(s) has a pronounced jump at e = e*. 

We have tested our approach in the Fee case and have got very good 
agreement with the computed values of D 2 from our linear fittings and the 
values reported in H . In figure 2b we show a typical curve used to estimate £*, 
and in figure 3 we compare the values of the exponents A and D 2 as calculated 
in figure 1 and also the matching values of D\ from equation (|8p for the Fee 
case. From now on we use this procedure to estimate the exponents D\ for 
the other substitution Hamiltonians (|l|) considered here. 

Its now time to discuss our numerical results and details of their imple- 
mentations. The return probability C(T) was calculated by direct diagonal- 
ization of the Hamiltonian and eigenfunction expansion of the initial state; 
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we used bases of size N m 1 x 10 3 and checked some results with bases of 
size N m 2 x 10 3 . The initial condition was always concentrated on the centre 
of the basis no and we have followed its time evolution until time Tf for which 
the modulus of the amplitude at one of the border sites reaches 1 x 10~ 6 . 
For the calculation of D 2 we have considered subdivisions of the spectrum 
in subintervals of size e ranging from the least eigenvalues spacing (we disre- 
garded multiple eigenvalues) up to 1CT 2 . We could seldom conceive a linear 
behaviour in such log — log plots in both D 2 and A cases in order to ex- 
tract faithful exponents, so that we were left with the task of finding e* and 
computing only D 2 . 

In figure 4 we present a summary of our main numerical results, i.e., the 
values of D 2 for some substitution potentials as function of the potential 
intensity A. Since we are not aware of any complete rigorous spectral classi- 
fication for Hamiltonian (|l|) in the cases of RS and PF sequences, we have 
also used Tf as indication of any possible (de) localization transition; this is 
the reason for the restriction of the RS case to A < 1.7; for all substitu- 
tion sequences we have found Tf < 10 3 for A < 2, but for RS Tf jumps 
from Tf pa 10 3 for A = 1.7 to Tf pa 10 5 for A = 1.8, which characterizes 
absence (at least numerically) of extended states. If non-localized states are 
present their "amounts" suffer a drastic reduction at A ~ 1.8 so that we have 
not detected them. Recall that in |25|, it is argued that all states of RS 
Hamiltonian should be localized for A > 2.0. 

Since no such sharp transition in Tf was found for the PF Hamiltonian, 
its values of D 2 are close to the corresponding values for PD, and both 
substitution sequences have point autocorrelation measures, we conjecture 
that the PF Hamiltonian (p]) has singular continuous spectrum for ^ |A| < 
2 (maybe also for any A ^ 0) as PD Hamiltonian does [18|. 



Besides the above conjecture we see from figure 4 that for all sequences 
the exponent D 2 decreases as A increases (as physically expected). Since 
different exponents were found despite the proven singular continuous spectra 
of Fee, TM and PD Hamiltonians, we see that D\ is able to discern these 
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operators. The values of for the RS are very close to the TM case, 
but in principle one would not expect this since the autocorrelation measure 
of the RS sequence is Lebesgue measure (the same for random sequences). 
Notice that the spectral classification of the underlying sequence generating 
the potential does not reflect exactly in D 2 , since these values for PF and PD 
are below the corresponding ones for TM and RS, the latter sequences being 
considered "more random" than the former ones. Notice, however, that for 
the most popular sequences Fee, TM and RS we have found agreements 
with the classification through the sequence spectral type, i.e., this order 
implying decreasing values of D 2 . 

The higher values of the exponents for RS compared to PD indicates 
that the presence of extended states mixed with numerically found localized 



states p5| , [26| does not necessarily imply lower exponents D 2 . 

Summing up, we have found that for substitution potentials in general the 
dynamical exponent A and the correlation dimension D 2 are difficult to be 
obtained from direct linear fittings (at least with the basis sizes we used; we 
suspect this is a consequence of the quantum intermitency and multiscaling 
in time of the dynamics fl2, |i3|Q and we have proposed D 2 as an effective 
exponent, which has recovered D 2 in the cases it can be directly obtained. 
We then computed D 2 for some substitution potentials and have not found 
a complete qualitative agreement with the spectral type of the substitution 
sequences themselves, i.e., more random the sequence smaller D 2 . Only for 
RS we have got indications of a spectral transition from extended (critical) 
to localized states, although its values of D 2 are higher than those for PD 
and PF cases. 
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Figure Captions 



Figure 1 

a) Log- log (base 10) of the return probability (dashed line) C(t) for the 
FCC potential with A = 1.0 versus time. The slope of the straight line fitting 
(full line) corresponding to A is indicated. 

(b) Log-log (base 10) of 7* (dashed line) for the Fee potential with 
A = 1.0 versus e. The slope of the straight line fitting (full line) corresponding 
to D 2 is indicated. 

Figure 2 

Same as in figure 1 but for the PD potential with A = 1.6. The arrow in b) 
indicates e*; the first point at left in b) corresponds to the least eigenvalue 
spacing, for which 7* = J2k \ a k\ A - No linear fitting is shown. 

Figure 3 

Scaling exponents D 2) D\ and A for various potential intensities A of the 
Fee potential. 

Figure 4 

Effective scaling exponents D\ as function of A for some substitution 
potentials. 
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